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We developed a microscopic approach to calculate the sample-to-sample fluctuations of trans- 
mission distribution in disordered conductors. This bridges between Green's function and random 
matrix theories of quantum transport. The results obtained show that the correlations of trans- 
mission eigenvalues obey universal Dyson statistics at small separations between eigenvalues being 
non-universal otherwise. The results facilitate an easy computation of different physical quantities 
and impose important constrains on a future imaginable theory of quantum transport. 



The concept of universal conductance fluctuations [1] proves to be one of the most interesting concepts in mesoscopic 
physics. Quantum interference results in conductance of a disordered coherent conductor fluctuates from sample to 
sample by a value of the order of Gq = e 2 /(2iTh). The exhaustive qualitative description of the phenomenon has been 
elaborated on the basis of Green's function technique. [2] 

The alternative explanation of UCF has been proposed by Imry [3], who asserted universal explanation of the 
phenomenon. The conductance of a mesoscopic sample can be expressed as a sum of the eigenvalues of the transmission 
matrix square, T = t^t, 
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!^ Basing on works by Pichard, [4] Imry showed that the assumption of Wigner-Dyson statistics of transmission 
^eigenvalues yields the correct order of magnitude of UCF. 

The pioneering paper by Imry has triggered persistent attempts to formulate the theory of quantum transport as a 
1 ~^kind of random matrix theory (RMT). This theory would produce the quantitative description of transmission matrix 
being based on a few general postulates. [5] However such a theory was doomed to fail. Green's function methods 
clearly showed that UCF are not precisely universal. Their magnitude depends on the shape of the sample and its 
^^dimensionality. Universal RMT would never have such input parameters. 

^ The puzzle is that the RMT can be made working. It may correctly describe quasi-one-dimensional geometry [6,7] 
QQand chaotic cavities. [8] Why? How could it be that the same phenomenon combines universal and non-universal 
Cofeatures? Are there ways to an improved RMT which would able to describe a concrete physical situation? 
IT) Below we answer these questions checking the Imry conjecture by the recently formulated [9] Green's function 
^^method. The short answer can be formulated as follows: correlations between transmission eigenvalues obey Wigner- 
-Dyson statistics provided their separation is much smaller than unity. At separations of the order of unity their 
^correlations become dependent on a concrete geometry of a sample. When calculating the fluctuation of conductance, 
^one sums up over small and moderate separations. This is why the final answer combines universal and sample- 
■^Jdependent features. The effective interaction between transmissions that corresponds to the correlations calculated, 
Cnwould be an input parameter of the imaginable improved RMT. 

O One can draw a parallel between statistics of transmission eigenvalues through an open systems and statistics 
of energy levels in a closed system. The latter was shown [10,11] to be universal at small energy separation and 
sample-dependent at energy differences = E t h- 

To proceed, we make use of multicomponet Green's function method proposed in [9] .It expresses generating function 
of momenta of the transmission distribution 

= Tr[ ?-——) (2) 

V ; L l - sin 2 ($/2)T ; V ; 

in terms of an artificially constructed Green's function which is a 2 x 2 matrix, energy. The Green's function introduced 
can be readily averaged over impurity configurations with using the conventional technique. [12] 

To find the correlations between transmissions, we have to compute the average product of two Green's functions 
of the kind. The technical details will be reported elsewhere. [13] The corresponding diagrams preserve the common 
cooperon diffusion structure [1]. However, those are difficult to handle due to the extra matrix structure and coordinate 
dependence of the average Green's function. The compact answer is possible to obtain for the expressions of the form 
(2). Using a common diagrammatic trick outlined in [14], we obtain 
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sin$sin$' <F($)F(&) >= 4— — (In det(l - K d ) + In det(l - K c )). (3) 

The operator K d ^ is actually an integral kernel corresponding to a diffusion (cooperon) ladder section. For a 
common case of contact impurity potential characterized by the elastic scattering time r it reads 

K d aliCd (x, x'; $, $') = -L < G °c( X} x >. $) >< G ad {X) x >. $ /) >; (4) 

K btCd (x, x'; $, $') = -L < G ac (x, $) >< G da (x', $') >, (5) 
t(x) 

with indexes labeling spin and the additional matrix structure 

To solve the problem, we have to find the eigenvalues of the operators 1 — K d ( c \p. It is possible to derive diffusion- 
like equations for those. These equations can be solved for a slab geometry of the conductor. In this case p can be 
readily expressed in terms of common diffusion and cooperon mode eigenvalues. 

These modes have been intensively studied, [12], and I shortly remind the reader the results. Let the slab dimensions 
be L x ,L y ,L z , x being transport direction. First, the discrete diffusion and cooperon modes can be labeled by three 
integer numbers n xyz , n y z ranging from to oo and n x starting from 1. The eigenvalues normalized by Thouless 
energy D/L 2 are given by ( n = (n x , n y , n z )) 

en = ^(nl + (^r + (^f). (6) 

ljy ljy 

It is easy to take into account spin effects and dephasing by magnetic field. To describe those, we introduce the 
following dimensionless parameters: r/so = Dtso I L x , Tso being characteristic time of spin-flip process, and r/ = 
eH L x L y I cTi, H being magnetic field || z. We can take different rj for different thus describing parametric correlations 
of Then all diffusion modes acquire an extra shift (ry — r/) 2 while the cooperon ones are shifted by (ry + r/') 2 . If 
spin-orbit interaction is taken into account, for each branch there are three degenerate eigenvalues corresponding to 
spin triplet and the one corresponding to spin singlet. The triplet eigenvalues acquire the same extra shift for all 
branches, e*">' e * = £ sin s let + Vs o- 

Due to extra matrix structure of the Green's functions in use there are two eigenvalues /j, corresponding to each 
either cooperon or diffusion mode, 



± _ \ 2 _ ($±$Q 2 



where mean free path / = s/Dt, m labeling all abovementioned modes. Substituting this to 3, we find 

sin* sin*' < F^)F(^) >= 4—— £ [ln( £m - 1 \ ) ) + ln(e m - 1 \ ' )]■ (8) 

modes 

Finally we express the answer in terms of density-density correlator. It is convenient to consider the density in 
new coordinate A related to transmissions by T = l/sinh 2 A. Using identity 2irp(X) = F(2X + it) — F(2X — it) that 
immediately follows from 2, we obtain 

< p(X,i 1 )p(X',i 1 ') >= -L V Z( v ^+7r) + Z( v ^-7r)-2Z( v ^); (9) 

Z7T Z z — 4 

modes 

Z(x; X, A') = M (x + i( l +x , ))2 + (x + i ^_ xl)) 2 )- 

This relation provides an exact microscopic answer for transmission correlations and is the main result of the present 
work. It describes the effects of dimensionality and crossovers between different statistical ensembles. 
Below we analyze the physically meaningful consequences of 10. 

First let us consider correlations at small separations between A's, |A' — A| <C 1- Them the sum is dominated by 
the second term corresponding n x = l,n y = n z = 0. Manipulating r),r)so allows to consider limits of different pure 
statistical ensembles and to reconcile that 

< ^ A, ) >= -^ Re (A-A' + iO) 2 ' ( 10 ) 
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s being an extra degeneracy of transmission eigenvalues not lifted by the random scattering. This coincides with a 
known result of Dyson [15] provided |A — A'| exceeds the average spacing between transmissions, sGq/G, and thus 
proves Imry conjecture. 

The present theory is not able to give a detailed form of the correlator at separations comparable with the average 
spacing, albeit we believe it coincides with Dyson expressions. This is a common drawback of theories perturbative 
in G/G q , like, for instance, [10]. At small separations, the fluctuations of the mode n x = 1, n yz = become big, like 
the fluctuations of the mode n xyz for an isolated system. [10]. For an isolated system, the answer has been gotten 
by non-perturbative supersymmetry technique. [11] The development of supersymmetry technique for diagrammatic 
approach [9] that has been started [16] recently, will do the same for transmission. 

If r) ^ rj' <C 7r the contribution of divergent mode 

< p(X, rj)p{\' ,r}') >= -^ (A _ A , + i( J ? _ J? , )/27r)2 > (11) 

reconciles universal parametric correlations [17] that correspond to Brownian motion of transmissions. [18] 

Thus we have shown that the transmission statistics at small separations obeys that one of Wigner-Dyson ensemble. 
The rest of terms in 10, however, do not have universal form, that explains why the fluctuations of the quantities 
integrated over transmission distribution, like conductance, do not obey universal relations. 

The question of interest is if one can improve RMT in such a way that it accounts for details of the sample geometry 
and describes crossover between ensembles. The answer seems to be positive. 

The proven way of such an improvement is so-called DMKP method [19] that seem to work for quasi-one-dimensional 
systems. [7] Indeed, for quasi-one-dimensional case, and for pure ensembles, ^JT^ = irn, sum in 10 become telescopic 
and can be readily done. The result coincides with the one obtained in [6] by DMKP method, 

< 'W'M >= ^ Re[ (A-A'W " (A-aW ' (12) 




FIG. 1. The Fourier component of effective eigenvalue interaction for ( from the uppermost to lowermost curve ) wire, 
Wigner-Dyson interaction, square, and cube. 

To become more optimistic, let us note that the correlator 10 possesses the following property, 

< p(X)p(X') >= K(\ - A') + A'(A + A'). (13) 

In the Green's function framework, it is impossible to find a physical reason for this. RMT gives a bright pattern to 
comprehend such a symmetry: consider a system of classical particles confined in the region < A. The correlator of 
such form emerges if one assumes pair interaction U(X — A') between the particles and their "image charges" mirroring 
the particles at negative A. [6] There is a simple relation between Fourier components of the interaction and K: 

u ^ = m) (14) 

In Fig.l I have plotted the Fourier component of the interaction for a pure ensemble and for different geometries 
(wire, L x L zy ; square, L x = L y L z ; and cube, L x = L y = L z ) along with Wigner-Dyson interaction. 
At moderately large k (small separations) for all geometries interaction rapidly converges to Wigner-Dyson values, 
whereas its small k behavior drastically differs. For a wire, U(k) — ► k~ 2 , and the eigenvalue repulsion exceeds the 
Wigner-Dyson one. For 2D geometry, U(k) — ► 2L X / L y \k\. For 3D case, U(k) — ► irL 2 /L y L z + c\k\. This corresponds 
to the following asymptotes at big A: 
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U(\)->-\\\ for ID 
U(X) —> —2L x /itLy In |A| for 2D 
U(\) -+ 2(L X /L y + L x /L z - 1)/A 2 for 3D. 



(15) 
(16) 
(17) 
(18) 



It is interesting to note that in 3D the asymptotical interaction can be attractive for sufficiently wide slabs. 

In conclusion, we have microscopically checked Imry conjecture for transmission matrix explicitly showing where 
and why Wigner-Dyson approach does not work. We have calculated the interaction that shall be postulated in an 
imaginable RMT that describes realistic diffusive conductors. 

I am indebted to C. W. J. Beenakker for his persistent interest to this work, teaching me RMT and many instructive 
suggestions made. I appreciate interesting discussions of the results with G. E. W. Bauer, K. Frahm, and B. Rejaei. 
This work is a part of the research program of the "Stichting voor Fundamenteel Onderzoek der Materie" (FOM), and 
I acknowledge the financial support from the "Nederlandse Organisatie voor Wetenschappelijk Onderzoek" (NWO). 
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